Eliminating nonlinear phase mismatch in resonantly enhanced 4- wave mixing 
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Resonantly enhanced four wave mixing in double-A systems is limited by ac-Stark induced non- 
linear phase shifts. With counter-propagating pump fields the intensity-phase coupling has minimal 
impact on the dynamics, but it is of critical importance for co-propagation. The nonlinear phase 
terms lead to an increase of the conversion length linearly proportional to the inverse seed inten- 
sity, while without nonlinear phase-mismatch the scaling is only logarithmic. We here show that 
the ac-Stark contributions can be eliminated while retaining the four-wave mixing contribution by 
choosing a suitable five level system with appropriate detunings. 

PACS numbers: 42.50.Gy, 32.80.Qk, 42.50.Hz 



Ever since the cancellation of resonant linear absorp- 
tion and refraction via electromagnetically induced trans- 
parency (EIT) jjj was first demonstrated, quantum and 
nonlinear optics have successfully been exploring the con- 
sequences. Many interesting effects have been proposed 
and investigated Q . One important application of EIT 
is optical frequency mixing close to atomic resonances 
where it allows making use of the resonantly enhanced 
nonlinear interaction without suffering from linear ab- 
sorption and refraction. It has been predicted that EIT 
could even lead to a new regime of nonlinear optics on 
the level of few light quanta || ||, |j| . 

In this paper we consider one particular EIT-based 
scheme, namely resonantly enhanced four wave mixing 
in a double lambda system as shown in Fig. [I]. The two 
fields with (complex) Rabi frequencies f^i and Q 2 are ini- 
tially excited and form the pump fields, while the other 
fields with (complex) Rabi frequencies E\ and E 2 are 
generated during the interaction process. and E\ are 
taken to be exactly on resonance while the other two are 
assumed to be detuned by an amount A. A finite detun- 
ing A, large compared to the Rabi frequencies, Doppler 
broadening and decay rates from the excited states, is 
necessary to maximize the ratio of nonlinear gain to lin- 
ear absorption. Decay from the two lower levels is con- 
sidered to be negligible. Because of energy conservation 
all fields are in four-photon resonance. It can be shown 
furthermore that the contributions of the resonant tran- 
sitions to the linear refractive index vanish if the fields 
are pairwise in two-photon resonance. Phase matching 
will thus favor two-photon resonance and we assume that 
this condition is fulfilled. Resonant four- wave mixing 
has been analyzed both theoretically and experimentally 
with co-propagating as well as counter-propagating fields 

1,0,1,1,0,0, 11 1|. 

Associated with the finite detuning A are ac-Stark 
shifts which lead to intensity dependent dynamical phase 
shifts of the fields. These phase shifts are of minor 
consequence in the case where the fields are counter- 
propagating pJ. They do have a detrimental influence, 
however, for co-propagation. In the following we will con- 
centrate on the latter situation and show how to elimi- 
nate these terms leading to a considerable improvement 



FIG. 1: Parametric amplification in a generic double-A sys- 
tem. 



of nonlinear frequency conversion. 

The standard method to describe the wave mixing 
process in a resonant medium is to derive density ma- 
trix equations for the atomic system, solve them in the 
steady-state, i.e. assuming adiabaticity, and insert the 
resulting expressions into the Maxwell-Bloch equations. 
This yields four equations of motion for the fields in the 
slowly varying amplitude and phase approximation. The 
field equations can then be further broken down into a 
set of five coupled equations consisting of four amplitude 
equations plus the equation of motion governing the rel- 
ative phase between the fields. 

This procedure can be rather cumbersome, particu- 
larly when several atomic levels are involved. A much 
simpler way to derive the field equations is given by the 
Hamiltonian approach introduced in JlJ], which we will 
use in the following. This method makes use of the fact 
that the polarization P of the medium can be expressed 
as a partial derivative of the time-averaged interaction 
energy per atom H with respect to the electric field or, 
equivalently, the Rabi frequencies E^ 
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A similar expression holds for the drive field polariza- 
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tions with E replaced by O. Here (...) denotes quantum- 
mechanical averaging, rfj the dipole matrix elements of 
the corresponding transitions, and N is the atomic den- 
sity. Hence introducing moving coordinates (C,t) with 
£ = z — ct one can directly obtain the field equations 
of motion in the slowly-varying amplitude and phase ap- 
proximation: 
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where r/^ = Ndfuji/(2hceQ). The evaluation of the right 
hand side of (|^) is particularly simple if an open-system 
model can be used to incorporate decay in a complex 
Hamiltonian H and if the atoms adiabatically follow the 
dynamics of the fields. If the atoms are initially in an 
eigenstate of H with eigenvalue A then (if) can simply be 
replaced by A as the two are equivalent. Thus knowledge 
of the eigenvalues of the single-atom interaction Hamil- 
tonian is sufficient to directly derive the field equations 
of motion. 

In the basis (|1) |2) |3) |4)) T the system shown in Fig- 
ure [l] can be described by the complex interaction Hamil- 
tonian 



H = h 



Taking O as a characteristic magnitude of the Rabi fre- 
quencies involved, to second order in O/A the relevant 
eigenvalue of (W is given by 
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This eigenvalue corresponds to the state asymptotically 
connected to |1) at t — > — oo. That is, the eigenstate 
associated with Ao corresponds to the ground state |1) 
for vanishing E\ and E 2 - If the pump fields change suf- 
ficiently slowly we may assume that all atoms stay at all 
times in this eigenstate and (H) in @ can be replaced 
by Aq • This yields the following equations of motion pit] 
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where k = 3NA 2 j/8ir, with 7 being the decay rate from 
the upper levels, where for simplicity we have assumed 
7i ~ 72 = 7- Note that there are no linear absorption 
terms in (||) - @, despite the presence of the decay terms 
7. Thus, the process is a parametric one, and the total 
energy of the electromagnetic fields is conserved. One 
furthermore finds that the equations have the following 
three constants of motion: 
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This allows the problem to be reduced to two degrees 
of freedom, one corresponding to the exchange energy 
between the fields and the other to the relative phase 
<p = <f>m + <f>m - 4>ei — 4*E2 between the fields. 

The terms in the second line of each equation in (|| — 1|) 
are ac-Stark induced, intensity dependent phase terms. 
They have a considerable impact on the dynamics, partic- 
ularly in terms of the conversion length, i.e. the distance 
required for one of the pump modes to attain maximum 
transfer into one of the generated modes. To see this we 
solve (||) — (|J) analytically, using the constants of mo- 
tion and the methods described in |l5|, [l6|. We will not 
give details of the derivation nor the full solution here, 
however, as they are not very instructive for the present 
purposes. 

We consider the case where the two pump fields are 
initially of equal intensity, as are the two seed fields, so 
that Ei = E 2 = E and Ox = 2 = 0. If we define the 
seed field parameter by e = |_E(0)| 2 /|O(0)| 2 and denote 
the initial relative phase difference between the fields as 
tpo then, in the limit e< 1, the conversion length L and 



conversion efficiency e = (| E n 
given by: 
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One immediately notices that in order to obtain full con- 
version as e — > it is necessary to carefully choose the 
initial phase ipo- Secondly the conversion length scales 
as L ~ l/\/e, that is, inversely in the seed field ampli- 
tude £7(0). Consequently for small, let alone vacuum, 
seed fields, the conversion distance will rapidly become 
infeasibly long. 
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On the other hand, if the phase terms in (|| — 1|) were 
to be omitted, the quantity 

Rein^E^E?,) = in^E^E^ cosp, (14) 

is another constant of motion. In this case, if one of the 
generated fields E\, Ei vanishes initially, i.e. if only one 
of these fields is seeded, then this constant of motion must 
be zero. This indicates that the relative phase <p can only 
jump discontinuously between ±7r/2, which occurs only 
at the end of each conversion cycle, when at least one 
of the field amplitudes vanishes. Thus in this case the 
phase is essentially decoupled from the evolution of the 
field amplitudes. This makes a considerable difference to 
the dynamics. From an analytic solution one finds for 
the conversion efficiency e and the conversion length L 
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where we have taken < cos <po < 1 for simplicity. Thus 
if the phase terms are not present, it is possible to obtain 
full conversion in the small seed field limit, regardless 
of the initial phase condition. Secondly, the conversion 
distance scales only as — loge, and the situation is com- 
pletely different to the previous case — the conversion 
length will always remain short. 

As an illustration of conversion distance dependence on 
seed field intensity we have calculated numerical solutions 
to Eqs. (||) — ((§]) with and without the phase terms, and 
without making the approximation that e<l. We have 
assumed Ey = E%, f2i = f2 2 and tpo — ir/4. The results 
are shown in Figure ||. 
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FIG. 2: Initial field intensities taken as |^i(0)| 2 = l^2(0)| 2 , 

|J*i(0)| 2 - " 
e=\E! 
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with the seed field intensity defined as 



Thus, given the obvious advantages inherent in the 
omission of these phase terms, the question naturally 
arises: does there exist a situation in which these terms 
can be made to vanish? 



It has been shown by Harris pTij that in a system of 
parallel A transitions with different excited-state energies 
there exists an optimum detuning such that the non- 
linear index of refraction vanishes. A similar idea can 
be applied here. Noting that both parts of (Ji|) are lin- 
ear in A, but only the first part is linear in each of the 
fields, suggests a method for canceling the phase terms 
and at the same time retaining the nonlinear interaction 
part. To see this, consider the five level system shown 
in Figure |^. Here we have introduced hyperfine angu- 



Pi(F=l) 




FIG. 3: Modified double-A system. F and ttif represent the 
total and z-projection angular momentum for the atom, in- 
cluding hyper-fine interaction. 



lar momentum terms, level |3) in the original scheme 
has been split into two and the pump beam 0,2 tuned 
to a point midway between them. Within the approx- 
imation that the electronic parts of the wave functions 
describing the two P\/ 2 states are identical, as are the 
electronic parts of the two S1/2 states, we find that the 
coupling strengths for the |5i/ 2 ) — > I-P1/2) transitions 
all have the same magnitude. Crucially, however, the 
\S 1/2 , F = l,m F = -1) -> |P 1/2 , F = 1, m F = 0) tran- 
sition has opposite sign to the other three. The Hamil- 
tonian thus becomes 
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We now find that to second order in f2/A the lowest 
eigenvalue of (p7|) is 



A„ = I 



\E!\ 



(18) 



We see that the phase terms responsible for the increase 
in the conversion length and sensitive dependence of the 
conversion efficiency on the initial phase are indeed ab- 
sent. It should be noted that under more general condi- 
tions, e.g. if the dipole moments to the F — 1 and F = 
manifolds are not equal, other values of the detunings 
need to be chosen. 
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In summary we have shown that the nonlinear phase 
contributions arising in resonant forward four-wave mix- 
ing due to the ac-Stark effect can be exactly eliminated if 
a five state system with appropriate couplings and detun- 
ings is used. We have derived a simple effective Hamilto- 



nian for this system and shown that the conversion length 
scales only logarithmically with the inverse seed intensity 
whereas with the phase terms present conversion length 
scales linearly. 
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